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Abstract 

Background: In our Institute PET/CT images are used for detecting, staging and monitoring various 
malignant tumours. Standard Uptake Value (SUV) is now common place in our clinical PET/CT 
oncology. Moreover, PET images are used for target volume definition in radiotherapy applications. 
Even if PET represents a powerful diagnostics tool, the quantitative data extracted from PET are 
affected by the limited resolution of the system. The aims of this work is to analysed by a phantom 
study the accuracy of the data quantified by our tomograph and the validation of the segmentation 
method used for radiotherapy applications. 
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Methods: The phantom study was performed on a PET/CT system using the IEC Phantom with 
spheres filled with different sphere/background activity ratios (Rs/b). Sphere volumes were determined 
using an adaptive thresholding method. SUV assessment was evaluated measuring the SUV values in 
each sphere and the percentage error respect to the true value. Hot Contrast Recovery Coefficients 
(HCRC) were also measured. 

Results: A linear relationship between thresholds and volumes was observed for volumes up to 10 mL. 
Between 10 mL and 5.5 mL thresholds decrease reaching a minimum at 1.1 mL. Lor volumes <1.1 mL 
thresholds increase exponentially. No dependence on the acquisition time was observed but thresholds 
depend on sphere volumes, Rs/b and smoothing filter. SUV is quantified for volumes up to 2.5 mL 
with an underestimation of 10%, for smaller dimensions SUV values are underestimated up to 80%. 
Hot Contrast Recovery Coefficients range from 87% to 16%. 

Conclusions: Lor objects with very small volume (<2.5 mL) the SUV values are affected by a 
significant error (up to 80%). Lrom the clinical point of view this means that very small lesions with 
low measured SUVmax might represent a false negative. Moreover, the limited PET resolution 
influences lesion segmentation: the adaptive thresholding method is an useful tool for tumour 
boundary definition but it may provide unreliable results for volumes less than 2.5 mL. 

Keywords: Positron emission tomography; Quantitative approaches; Segmentation; SUV 


1. Introduction 

PET/CT has revolutionized medical diagnosis in many fields, by adding precision of anatomic 
localization to functional imaging. The benefits of PET/CT images in the entire spectrum of cancer 
care, from diagnosis and treatment planning to evaluation of treatment responses, have come into focus 
in recent years as the technology has become a standard practice [1,2]. Moreover PET images are 
increasingly used in Radiotherapy to define the tumour boundaries [3]. High precision Radiotherapy 
requires accurate tumour volume definition. Target volume contouring is usually based on the 
anatomic information acquired with Computerized Tomography (CT) or Magnetic Resonance (MR). 

In the last years the introduction of PET in combination with CT provides a better definition of the 
Biological Target Volume (BTV) and therefore improves tumour localization and delineation [4,5,6]. 
However partial volume effects [7] caused by the limited spatial resolution of the imaging system 
affects quantification of PET data. These effects and also the parameters used in image reconstruction 
influence the segmentation methods used to define the tumour boundaries and the quantification of the 
activity uptake in terms of Standardized Uptake Values and as consequence the tumour assessment. 
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For all these reasons, before using PET images in Diagnostics and Radiotherapy applications, it is 
important to analyze the influence of partial volume effects and image reconstruction parameters on 
PET data quantification. 

In our Institute PET/CT images are used for detection and staging of various malignant tumours, 
monitoring of their response to therapy and also for Radiotherapy tumour volume definition. For 
detection, staging and monitoring SUVmax values are used to define the pathologic relevance of a 
lesion. 

For Radiotherapy applications the definition of the biological target volume (BTV) relies on automatic 
segmentation of PET images using an adaptive thresholding method based on tumour/background 
ratios [8,9,10]. 

Aim of this study is to define with a phantom study the actual contouring thresholds to be used in 
tumour boundary definition and to evaluate the accuracy of PET data quantification in our tomograph. 

2. Materials and Methods 

The study was performed on a DISCOVERY TM 710 PET/CT scanner (Ge Healthcare) [11]. The 
system has lutetium orthosilicate detectors covering an axial FOV of 16.2 cm and a transaxial FOV of 
70 cm in diameter. PET data were acquired in 3D-mode. The image matrix was 256 x 256 with 1.56 
mm pixel. The PET image slice thickness was 3.75 mm. The CT scans were performed with 120 kV, 
400 mA, 0.5 sec tube rotation, and 3.75 mm slice thickness. CT data were used for attenuation 
correction. The acquired data were corrected for decay and dead time. 

Measurements were performed on a cylindrical phantom which is a modified version of the IEC Image 
Quality Phantom [12] (Figure 1) containing three more spheres for a total of 9 spheres with volumes 
ranging from 0.14 mL (diameter of 0.6 cm) to 98 mL (diameter of 5.7 cm). 



Figure 1 . Sections of the IEC Image Quality Phantom. 
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The phantom simulates a human thorax with lesions of different sizes and activity concentrations. The 
sphere/background activity ratios (Rs/b) ranged from 2 to 73 for a total of 11 phantom sets. For each 
phantom set PET scans in 3-D list mode were performed with different acquisition times (2, 3, 4, 5 
minutes) at the aim to evaluate the influence of the acquisition time on data quantification. 

Images were reconstructed with the following parameters [13,14]: 

• OSEM algorithms (Ordered-Subset Expectation Maximization) including TOF (Time of Flight) 
and PSF (Point Spread Function); 

• 3 iterations and 18 subsets; 

• a z-filter with a full width at half maximum (FWHM) of 4 mm. 

Each acquisition was elaborated with three transaxial filters of 4 mm, 6 mm and 8 mm (FWHM) at the 
aim to evaluate the influence of the filter on data quantification. 

A total of 132 datasets were obtained and analysed. Each dataset is identified by the Rs/b value, the 
acquisition time and the z-filter. 

Before the phantom experiment a cross calibration between the PET camera and the dose calibrator 
used to measure the 18 F activity was performed with the aim of minimizing the differences in activity 
measurement between the two systems [15, 16]. 

For each dataset the following physical figures of merit were considered: 

• Contouring thresholds to define the actual sphere volumes 

• Maximum and medium 18 F concentration in each sphere 

• SUVmax in the spheres 

• Hot Contrast Recovery Coefficients (HCRC) 

• Signal to Noise Ratio (SNR) 

• Contrast 

• Coefficient of Variance (CV) in the background 

The contouring thresholds were determined using a segmentation tool (PETVCAR GE Healthcare) 
based on an adaptive thresholding method where thresholds are defined as a percentage of the 
maximum activity in the sphere [17]. The sphere volumes were determined applying the threshold 
value that minimizes the difference between the actual volume and the calculated volume. A linear 
regression model was used to determine the relationship between contouring thresholds, Rs/b, sphere 
volumes and transaxial filter. 
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The SUV values were determined by the following formula: 


SUV= 


Act lol (kBq/ mL) 
Ac tadmm (MBq)/BW(Kg ) 


where Act vo , is the activity concentration measured in the volume of interest; Act admin is the total 
activity in the phantom corrected for the physical decay of 18 F and BW is the phantom weight (about 
10 kg). 

Definition of the Hot Contrast Recovery Coefficients is determinant to evaluate the activity 
underestimation in small objects due to partial volume effects [18]. The coefficients can be calculated 
using the following formula [19]: 

Measu red m eansphere activity— Measu red bkg activity 

ti C/v C = 

Known sphere activity— Knownbkg activity 


SNR, Contrast and CV characterize the image quality and they are influenced by the choice of the 
image reconstruction parameters and in particular by the value of the transaxial filter. 

The parameters were calculated taking into account the following relations [20]: 

SNR = [Act max - background] /SD 
Contrast = Act max /background (mean value) 

CV= SD/background ( mean value) (%) 

where Actmax is the maximum activity concentration value measured in the volume of interest, the 
background activity is determined as the mean activity in 10 ROIs placed in the background and the 
SD is the mean SD of these 10 ROIs. 

3. Results and Discussion 


3.1. Study Results 

The relationship between the measured thresholds and the sphere volume is linear for volumes up to 10 
mL while for smaller volumes the two variables show an exponential dependence. This is due to partial 
volume effects which are more significant for small volumes [21,22]. The observed trend is similar for 
all 132 datasets acquired. Figure 2 shows the threshold curves determined with an acquisition of 3 
minutes and with a transaxial filter of 6 mm filter. 


The range of threshold values is between 30% -f- 87%, with the lowest values for high Rs/b (Rs/b = 73) 
and the highest values for low Rs/b (Rs/b = 2). Different threshold trends were observed for volumes 
higher than 5.5 mL and volumes lower than 5.5 mL. 
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For volumes up to 5.5 mL the threshold values are included in the range of 37% 4- 42% for Rs/b 
between 70 and 8. In case of lower Rs/b (< 8) threshold values increase up to 56%. For volumes lower 
than 5.5 mL the threshold trend changes completely: a minimum value is observed for volume of 1.1 
mL and Rs/b in the range 8 4 73; this minimum shifts towards volumes of 5.5 mL for Rs/b lower than 
8. In case of very small volumes (<1.1 mL) a strong increase in threshold values is observed for all 
Rs/b. 



Figure 2. Measured thresholds versus sphere volumes for different background/sphere ratio (acquisition for 3 minutes and 
6 mm transaxial filter). 

Figure 3 shows the threshold trends as a function of the transaxial filter (4 mm, 6 mm, 8 mm) with 
fixed acquisition time of 3 minutes. The threshold values depend on the filter used on the image 
reconstruction, for a 8 mm- filter the threshold values increase of about 10% respect to an image 
reconstruction using a 4 mm-filter. This dependence is similar for all Rs/b and acquisition times. 

No dependence on the acquisition time was observed for threshold values except for the case of 2 
minutes where the low acquisition statistics causes a decrease of about 2% in threshold values. 
Therefore a 3 minutes acquisition time is sufficient to acquire an image with enough statistics to 
evaluate the patient data correctly. 


ISSN: 2057-3782 (Online) 
http://dx.doi.org/! 0. 1 7229/jdit.20 15-0216-013 


23 


Journal of Diagnostic Imaging in Therapy. 2015; 2(1): 18-29 


Ferrando et al. 


Ratio Sphere/Background: 73 



Figure 3. Threshold trend for image reconstruction with transaxial filter of 4 mm, 6 mm, 8 mm and Rs/b 73. 


The trend for volumes >5.5 mL was fitted with a linear regression model to define the equation 
relating thresholds, sphere volumes, Rs/b and the transaxial filter. From the fit the relationship between 
these quantities can be expressed as follows: 

Threshold % = + B l ■ FWHM fUter + B 2 ■ Sphere volume + B 0 (1) 


with A = 47.2 + 5, Bt = 1.78 ± 0.2, B 2 = 0.014 ± 2.1- 10' 3 , B 0 =23 ± 3 

For volumes between 2.5 mL and 5.5 mL the same equation was used to define contouring thresholds 
with an estimated error of ± 7% on the threshold values. For volumes < 2.5 mL partial volume effects 
severely affect activity quantification as shown in literature [23] therefore PET image segmentation 
using thresholding methods is not recommended since it can produces unreliable results. 

Relation (1) can be used to delineate the tumour volumes (> 2.5 mL) in clinical cases taking into 
account the signal to background ratio of the lesions. 

SUV quantification was evaluated considering for each dataset the maximum SUV value (SUVmax) 
and the percentage error on SUVmax respect to the true value. Figure 4 shows the SUVmax values for 
each Rs/b. 
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For all ratio Rs/b, the SUV is quantified for volumes up to 2.5 mL with an error of about -10%. For 
smaller volumes the tomograph underestimates the SUV values with an error of about -35% for 0.5 mL 
up to -80% for very small volumes (0.14 mL). Since the maximum SUV value is a parameter used in 
clinical assessment to define the pathologic relevance of a lesion, for a correct quantification of SUV 
the Hot Contrast Recovery Coefficients to be applied to obtain the actual SUV values were measured. 

As it is known partial volume effects influence the activity quantification in small lesions and therefore 
the SUV values. This is related to the scanner resolution. Each PET scanner has an intrinsic point 
spread function which is a profile of a point source in air. The width of this PSF is a measure of the 
spatial resolution of the scanner. For objects greater than two times the spatial resolution partial 
volume effects are negligible. If the objects are smaller than two times the spatial resolution a 
significant fraction of counts spills out of the reconstructed object. 

This can result in an underestimation of activity concentrations and overestimation of the object size. 
To compensate for the decrease in measured activity the recovery coefficients for hot spheres has to be 
calculated as a function of the object size. Table 1. Shows the HCRC measured for all Rs/b. Recovery 
coefficients range from 87% to 16% as the diameter of the spheres decreased. 


Sphere volumes (mL) 

Ratio 73 

Ratio 51 

Ratio 25 

Ratio 14 

Ratio 8 

Ratio 4 

98 

74% 

77% 

75% 

80% 

80% 

87% 

25.5 

69% 

72% 

70% 

74% 

75% 

82% 

11.6 

63% 

68% 

66% 

69% 

71% 

78% 

5.5 

59% 

65% 

63% 

68% 

70% 

73% 

2.55 

51% 

59% 

57% 

61% 

62% 

69% 

1.1 

50% 

51% 

50% 

55% 

53% 

57% 

0.5 

40% 

43% 

40% 

42% 

43% 

46% 

0.3 

29% 

27% 

26% 

28% 

29% 

39% 

0.14 

24% 

17% 

16% 

20% 

23% 

38% 


Table 1 . Hot Contrast Recovery Coefficients (HCRC) for different ratio Rs/b- 

SNR, Contrast and CV analysis are shown in Table 2. The higher value of SNR is obtained with filter 
of 8 mm. The Contrast and CV decrease when increasing the filter width. 
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Ratio 73 

Ratio 51 

Ratio 25 

Ratio 14 

Ratio 8 

Ratio 4 

SNR 

4 mm 

629 

554 

457 

228 

160 

83 

6 mm 

784 

761 

504 

258 

175 

97 

8 mm 

953 

849 

568 

310 

223 

115 

CONTI 

RAST 

4 mm 

72.80 

49.90 

25.33 

15.20 

8.50 

4.64 

6 mm 

68.20 

47.70 

24.22 

14.00 

7.85 

4.33 

8 mm 

65.70 

46.80 

23.50 

13.49 

7.57 

4.10 

CONTI 

RAST VARIABILITY (CV) 

4 mm 

11.4% 

8.8% 

5.2% 

6.3% 

4.7% 

4.4% 

6 mm 

8.6% 

6.1% 

4.5% 

5.0% 

4.0% 

3.4% 

8 mm 

6.8% 

5.4% 

3.8% 

4.0% 

3.0% 

2.7% 


Table 2. SNR, Contrast and CV versus FWHM transaxial filter (4 mm, 6 mm, 8 mm). 


3.2. Discussion 

Since in our Institute PET images are used to identify the Radiotherapy target volume and SUV values 
are used to classify the malignant grade of a lesion, aim of this study was to investigate the 
performance of our PET scanner in terms of SUV quantification and validate a PET segmentation 
method based on adaptive threshold algorithms. 

The study was performed using a cylindrical phantom containing hot spheres in a warm background. 
Different tumour/background ratio (Rs/b) were used to simulate concentrations of 18 F similar to those 
observed in clinical cases. The Rs/b range between 2 to 73, each phantom was acquired for different 
times (2 min, 3 min, 4 min, 5 min), the data were then elaborated including TOF and PSF algorithms 
and different transaxial filters (4 mm, 6 mm, 8 mm) to analyse the influence of the filter on data 
quantification. For each dataset the sphere volumes were determined applying adaptive thresholds 
calculated as a percentage of the maximum activity concentration in the spheres. For all ratios Rs/b the 
measured thresholds are linear for volumes >5.5 mL (sphere diameter of 22 mm). 

For volumes between 2.5 mL and 5.5 mL the linear trend is no more valid; in this interval thresholds 
decrease reaching a minimum for sphere volumes of 1.1 mF (diameter 13 mm). For volumes <1.1 mF 
thresholds increase exponentially because of partial volume effects. Thresholds vary in the range 30% 
■f 87% with lower values for high Rs/b (73) and higher values for low Rs/b (2). 
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The measured thresholds present a dependence on the transaxial filter used in the image reconstruction, 
increasing the filter value it is necessary to increase the threshold contouring a fixed sphere volume. 
Thresholds do not depend on the acquisition time; 3 minutes acquisition time are enough to reach a 
good statistics for data analysis. The trend of the measured thresholds for volumes above 2.5 m can be 
described with a linear regression model depending on three parameters: the sphere volumes, the ratio 
Rs/b, the filters applied in the image reconstruction. 

For volumes smaller than 2.5 mL the adaptive thresholding method is no more applicable because of 
partial volume effects therefore PET segmentation of lesions with dimensions < 2.5 mL seems to be 
inadequate. 

SUV quantification was investigated measuring in each dataset the maximum values in the hot spheres 
and the percentage error on SUVmax respect to the true value. For all ratio Rs/b, the SUV is quantified 
for volumes up to 2.5 mL with an error of about -10% while for lower volumes the tomograph 
underestimates the SUV values with differences of -34% for 1.1 mL up to -80% for very small 
volumes (0.14 mL). This underestimation of SUV in small volumes is caused by partial volume 
effects. Since the SUVmax value is a parameter commonly used to identify the malignance of a lesion, 
a quantification of SUV for lesion smaller than 2.5 mL is inadequate without taking into account a 
compensation of the decreased measured activity in the lesion. 

Lor this reason the Hot Contrast Recovery Coefficients for hot spheres were measured. Lor each ratio 
Rs/b the underestimation of the mean activity concentration as a function of sphere volumes was 
calculated. Recovery coefficients range from 87% (for sphere of 0.14 mL) to 16% (for sphere of 98 
mL). Lor each dataset the SNR, CV and Contrast characterizing the PET image quality were also 
analyzed. These parameters were evaluated as a function of the applied transaxial filter. The higher 
value of SNR is obtained with a filter of 8 mm. The Contrast and CV decrease with increasing the 
filter width. Lrom our analysis a filter width of 6 mm seems to be a good compromise to obtain a 
higher SNR and Contrast and lower CV. 

4. Conclusions 

The present study highlights that PET data quantification in terms of SUV can be strongly influenced 
by the limited PET resolution. Lrom the clinical point of view this means that very small lesions (<2.5 
mL) with low measured SUVmax (< 2g/mL) might represent a false negative. Moreover the limited 
intrinsic PET resolution and the choice of the image reconstruction parameters influence lesion 
segmentation of PET images. The adaptive thresholding method based on Rs/b used in this study may 
be an useful tool for tumour boundary definition but, as demonstrated, it may provide unreliable results 
for volumes less than 2.5 mL. 
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